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Abstract- In this paper a hand vein pattern authentication 
based on the frequency domain patterns correlation is 
presented. Correlation matching based on the hand vein pattern 
frequency domain spectrum is analyzed. In order to evaluate the 
system testing performance, a dataset of 50 persons of different 
ages above 16 and of different gender, each has 10 images per 
person was acquired at different intervals, 5 images for left hand 
and 5 images for right hand. In verification testing analysis we 
used 3 images to represent the templates and 2 images to test. 
Each image is matched with the existing 3 templates. We found 
that our system can operate at 72.72 percent genuine acceptance 
rate and a 0.047 percent false acceptance rate, and the 
corresponding threshold is 91. The system's testing EER is 10.19 
percent at threshold 87. For the 300*300 templates, the optimal 
threshold is 91 with FAR of 0.046 and FRR of 22.56. For the 
testing phase, matching the 200 test patterns with the 300 
templates of 100 hands, the optimal threshold is 91 with FAR of 
0.136 and FRR of 13.5 with EER of 3 percent at threshold 89. 
The system efficiency at this threshold was found to be 99.73% . 
K^^uwds- Biometrics, Verification, Hand Veins, FFT, Statistical 
Performance. 

1. Introduction 

AuthentLcating the identity of an individual using vein 
patterns was investigabed b^ several authors 11-5]. The 
infrared region is of special advantage since the sldn tissue is 
relatively transparent and the blood absorbs infrared light 
well. Hence, the image contrast is higher than in the visible 
area Since the arrival of fairly low cost CCD cameras and 
conputer power, it seems straightforward to try to consider 
these technologies 16-7]. There are many research attenpts 
for the extraction, segmentation and tracing of subcutaneous 
peripheral venous patterns 18-11], its main aim is to make 
data reduction and noise suppression for good diagnostic 
purposes and for making some quantitative measurements 
like lengths and diameters for the extracted vessel segments. 
These techniques based on mathematical morphology and 
curvature (vdns direction) evaluation for the detection of 
vessel patterns in a noisy environment Researchers in hand 
vein biometrics 11-5] had a satisfactory result for either 
verification or identification purposes, regardless of the 
difference in datasets size, methods, or vdn similarities used 
The vein tree detection stage includes four consecutive sub 
stages, which are hand region segmentation (i.e. region of 
interest localization and backgreund elimination), smoothing 
and noise reduction, local thresholding for separating vdns, 
and postprocessing. In literature of hand vein biometrics, 
Tanaka et al 13] have used phase oriy correlation to matoh 
two hand vein images, however he was working directly on 
the gr^scale pattern of hand vdn without extraction of that 
unique pattern In this p^Der, we propose a method for 
correlation of hand vdn patterns using frequency domain 
spectral matohing. 



II. Data Acquisition And Processing 

A monochrome frame-grabber is used to c^Dture an image 
of the back of a hand for conputer processing. Images are 
c^Dtured using a 320W X 240H pixels of 8-bits per pixel. A 
sample image is shown in Frg.l for a male hand A dataset of 
50 persons of different ages above 16 and of different gender, 
each has 10 images per person was acquired at different 
intervals, 5 images for left hand and 5 images for right hand 
The data set is for normal persons who do not conplain from 
any diseases such as arthritis. 




. Acquired image 



A. Hand Vein Image Processing Stages 



-bits per pixel 



This is the second stage in the Hand Vdn Verification 
System (HWS), which covers the detection of vein 
structures from the acquired irfered image for the back of 
the hand The vein tree detection stage includes four steps, 
which are hand region segmentation (i.e. region of interest 
localization and background elimination), smoothing and 
noise reduction, local thresholding for separating vdns and 
finally the postprecessing. 

B. Hand Region Segmentation 

Image segmentation is one of the most important steps 
leading to the analysis of processed image data. Its main goal 
is to divide an image into parts that have a strong correlation 
with objects or areas of the real world contained in the 
image. Binarization is the case of segmenting the image into 
two levels; object (hand region) and background; the object 
segment which is the region of interest (ROl) in white and 
the background segment in black as shown in Fig. 2. The 
algorithm used in the segmentation sub stage is an iterative 
method used for calculating and selecting an optimal 
threshold, which is used to segment the image into two 
distinct parts; hand and background [13]. We used this 
resultant binary image to calculate the center of gravity 
(COG) for our ROl (hand region). Then we translated the 
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grayscale hand region to the center of the image after 
assigning the background area to zero value pixels. Thus we 
completely localized, separated and centered the hand region 
for subsequent processing steps. 

lUU 

Fig.2. Segmentation results; (a) Input gray scale image (b) Binary image and 
(c) Output image after ROl determination and centering 

C. Smoothing and Noise Reduction 

Two ^proaches could be used for noise filtering. First 
^proach is using Gaussian smoothing filter. The 
disadvantage of Gaussian filter is its non-edge preserving 
ability, since it bluis the image with equal weights; also edges 
of the veins are bluned and conplebely diffuse aftjer several 
smoothing iterations. The second ^preach is an edge- 
preserving technique like nonlinear diffusion 116-17]; in 
which the image gradient was used to weight the diffusion 
process. Fig. 3 show results of the smoothing filters used in 
this work where we used a rnedian filter of 5* 5 rnask in Older 
to remove the hand traces from the acquired image then we 
used the nonlinear diffusion filter based on edge weighted 
diffusion in order to smoothen the image while preserving the 
vein edges. The smoothing and noise cleaning sub stage 
effect is shown in Fig. 3 for 5 iterations of noriinear 
diffusion, while the edges are not affected 



conputing the threshold for binarizing the central pixel. The 
result is shown in Fig. 4. 




Fig. 3. Effect of smoothing sub stage on the image line profiles 

D.Hand Vein Pattern Segmentation 

Specifically, hand van segmentation is to divide a hand 
vein image into a foreground (veins in the back of the hand) 
and a background (non-vessel areas). Segmentation methods 
can be divided into four groups, which are threshold-based 
segmentation, edge based segmentation, region based 
segmentation and segmentation by matching. In this study, 
the first thresholding method is adopted since it is 
computationally che^ and fast Considering that we want to 
process and study veins oriy, ^obal thresholding (i.e. single 
threshold for the whole image) is not a good technique for 
this purpose. A belter ^preach is to calculate the average 
around each pixel of the image in an area of NxN neighbor 
pixels and to use average value as a threshold value 111]. The 
local threshold process separates the vein pattern from the 
background; hence the desired vein image is extracted 
Experimentally we have chosen a 31x31 mask size for 




Fig. 4. Processed image (left) and its local thresholded image (right) 

E. Hand Vein Pattern Postprocessing 

It is demonstrated from Fig. 4 that the resultant binary 
hand vein contains some noise and un-sharp edges. We 
experimentally ^plied 5x5 median filter for inproving and 
validating the output binary hand vein pattem and for 
reducing the effect of these unwanted defects 113]. We also 
converted the vdn pattem into white in a black background 
which in this case the entire image. The final pattem after the 
post precessina sub stage is shown in Ficr. 5. 




Fig. 5. Hand vein pattern before (left) and after postprocessing (right) 

F. F re quency D omain Analysis of Hand Vein Patterns 

The resulted image after preprocessing was transformed 
using FFT. Example of FFT spectral image for a gr^scale 
pattem is shown in Fig. 6. Most of spectral patterns are star 
like sh^Des. Star r^s differ in size, sh^Des, and orientation 
about central axes. Information extracted from frequency 
domain from the star patterns can be later used in hand veins 
' domain patterns type classification 





Fig. 6. G rayscale vein (left), Pattern (M iddle) and FFT spectrum (right) 

G. Matching of Hand Vein Patterns in Frequency Domain 

After image acquisition and hand vdn extraction sub 
stages, we have a binary image contains the segmented vein 
pattem of the back of the hand The Fourier transform for 
this binary image was calculated using the Fast Fourier 
transform (FFT) algorithm; the 2-D Discrete Fourier 
Transform (DFT) is demonstrated in equation (1). We 
calculated the Fourier Spedium for the resulted patterns 
using equation (2) 113]. 



iux/M+vyl N) 



M-lN-1 
x=0 y=0 

\F{u,v)\=[R'iu,v) + l'{u,v)f" 



(1) 
(2) 
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This is suitable for the next and the final sub stage, the 
matching of hand vein patterns in frequency domain. As it is 
expected, the input for the matching sub stage is two hand 
vein spectral images like the one in Fig. 6 (right), the 
matching output is Yes (the two images are for the same 
pattern) or No (the input images is not correlated). We 
performed a normalized cross correlation between the two 
spectral images to get a normalized spectral similarity ratio 
between the two spectral patterns. Since frequency domain 
spectrum is invariant to translations in spatial domain. We 
address only rotation in registration technique to find the 
maximum spectral correlation [14]. We have already 
constrained our data acquisition system with the hand 
attachment in order to prevent any rotation. One of the two 
spectral images is remained stationaiy while we apply 2D 
transformation (rotation) on the other image in order to align 
it with the first pattern (Registration) to find the maximum 
correlation percentage between two hand vein spectral 
images [22-24]. 



Correlatwn{6) = ForaU{\X |* |T |*100)/Min(|X|',|r|' 
MaxCorrelation = Max{Correlation{0)) 



(3) 

(4) 



Where |X| is an image contains FT spectrum of the first 
test hand vein pattern, |T| is the second image contains FT 
spectrum of the transformed template hand vein pattern, * is 
the ordinary multiplication mathematical operator, Min 
(|X p,rrf ) is the minimum value of the summation of the 
spectnjm square between the two patterns to be niatx^ 

However Fujitsu lesearchets have piesenbed a 
contaclies^touchless palm vein authentication which is more 
convenient and non hygienic for users we restricted our 
designed prototype for a contact system type in order to 
simplify our matching phase and derive a dependent 
statistical measures for this proposed prototype. However 
accounting for scale (Age or hand to camera dstance) in the 
registration algorithm is sinple to inplement in frequency 
domain, we did not account for the scale in our matching 
metiiod In real systems c^Dturing the hand template at equal 
montiiy intervals after correct authentication is suggested in 
order to track the changes in the hand with age. The 
matching (similarity) percentage is calculated cross 
correlation between fhe two images spectrum We calculated 
the matching ratio for each transformation step then we 
choose the maximum ratio as the final matching ratio 
between the two irput hand van spectral patterns. In our 
inplementation and for saving time of matching we made 
the rotation parameter (rotation) steps equals five degree in 
frequency domain, getting the maximum matching ratio on 
this grid, finally we made a fine search (tune) to find the 
overall maximum correlation ratio. The result of matching 
sub-stage is shown in Frg 7, in case of correct true match it 
is demonstrated that the resultant pattem is correlated to the 
irput images and it is shown that the matching ratio is 94.00 
% (same person). A case of correct mismatch is shown in 
Fig 8, it is shown that the matching ratio is small as 73.43 % 
(different persons). 




Fig.8. Example of correct mismatch between different persons 



III. Results 



A. Statistical Results 



The system was tested over a database collected using the 
designed system consisting of 50 persons of df f erert age and 
gender for each 5 left and 5 right images were acquired The 
person was asked to put his/her right hand on the hand 
attachment frame and the system operator c^Dtures the first 
image for the current person right hand veins. Then the 
person replaces his/her right hand with the left one for 
acquiring the first image for the left hand vein pattern This 
process was repeated until we acquire five images for the 
right hand and five images for the left hand in different 
scenes (5 ninutes interval between every acquired image) 
independert of each other, i.e. ten images for each person 
We will prove in our statistical analysis that the harid vdn 
pattem is unique to some level for each person and for each 
hand Thus we considered as if we have 100 persons of 
which 5 images are in the dataset since we found that left and 
right hand vein images are differeit In order to find the 
dissimilarity threshold in correlation ratio between the 100 
subjects we have chosen oriy the first image pattem for each 
of the 100 subjects and for a correlation ratio threshold that 
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exceed 91%, we achieved 100% classification (Distinct 
pattern). For evaluating (testing) our system for verification 
purpose, all possible comparisons are made between the 
whole data. We matched each image from our data set with 
all the 500 hand vein images in our data set and then we 
recorded the matching ratios. We constructed the correlation 
matrix for representing the matching result between each 
image and all other images. We used the thresholded 
correlation matrix for calculating the efficiency of our system 
at each threshold. Also, we deduced that the optimal 
threshold at which the maximum efficiency occurs in our 
system. We performed statistical analysis for selecting an 
optimal threshold to get the highest system performance, by 
testing the system over the whole database - exhausting the 
whole database. To evaluate the hand vein performance we 
used measures of performance, which include: false accept 
rate (FAR), false reject rate (FRR), and efficiency. Fig. 11 
shows how FAR(%), and FRR (%) change with different 
thresholds. Our aim is to select an optimal threshold. We 
want a single criterion, such that it takes into account the 
maximization of true events (GAR, GRR) and minimization 
of error events (FAR, FRR). 

B. Verification Testing Results 

To obtain the verificatLon accuracy of our systern, each of 
the images was matjched with all of the images in the 
database after reducing it to 300 patterns (3 images for each 
hand instead of 5). A matehing is noted as a conect matehing 
if two images are from the same hand The total number of 
matohings is 20000. The prDbability distributions for genuine 
and impostor are estimatod by the conect and incorrect 
matohings, respectively, and are shown in Fig. 9. Frg 10 
depicts the corresponding Receiver Operating Charactoristic 
(ROC) curve, for all possible operating points. From Fig 10, 
we can see that our systom can operate at a 86.5 percent 
genuine acceptance rate and a 0. 136 percent false acceptance 
rate, and the corresponding threshold is 91. The system's 
testing equal enur rate is 3 percert at threshold 89. In this 
verification testing analysis we used 3 images to represert the 
ternplabes and 2 images to test Each image is matohed with 
the existing 3 tonplabes. The same analysis for the 300^300 
as in section A was repeated Frg. 9-11 ^ows the results. For 
the 300^300 tonplates, the optimal threshold is 91 with FAR 
of 0.046 and FRR of 22.55. For the testing phase, matehing 
the 200 test with the 300 tonplabes of 100 hands, the optimal 
threshold is 91 with FAR of 0.136 and FRR of 13.5. The 
systom efficiency at this threshold was found to be 99.73%. 
Fig. 1 1 shows FAR and FRR versus threshold 
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Fig 10. Verification test results the receiver operator characteristic curve. 
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Fig 9. V erification test results for genuine and imposter distributions 



Fig 11. FAR & FRR for template and test analysis 

IV. Discussion And Conclusion 

The designed system was tested for verification purpose 
oriy over a database collected with the designed system. 
Dataset for 50 persons of different age and gender of which 
ten images per person were acquired (five for the right hand 
and five for the left) in different scenes at different intervals 
and are independent of each other, i.e. ten images for each 
person Verification performance statistical parnmetets were 
estimated for the oveiBll system such as: Genuine Accept 
Rate (Sensitivity), Genuine Rgect Rate (Specificity), False 
Accept Rate (FAR), False Rgect Rate (FRR), Efficiency and 
Recdver Operating Curve (ROC). System testing 
performance (overall efficiency) was found to be 99.73% at 
threshold (matehing ratio) equal 91. At this maximum 
efficiency FAR of 0.136 and FRR of 13.5 were reported 
However the difference in methods, datasets, and algorithms 
that were found in the hand biometric work of [1-5, 21], our 
performance results are conparable. 
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